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We introduce a new kinematic event variable rriTGea which 
can provide information relating to the mass scales of par- 
ticles pair-produced at hadronic and leptonic colliders. The 
variable is of particular use in events with a large number 
of particles in the final state when some of those particles 
are massive and not detected, such as may arise in R-parity- 
conserving supersymmetry. 

1 Introduction 

When LHC experiments like ATLAS [HE] and CMS [3] begin to look for signs 
of supersymmetry (or any other models with large numbers of new particles), 
one of the first things they will want to probe is the mass scale associated 
with these new particles, if there are any. Single particles produced in narrow 
s-channel resonances will in general be easy to spot (one would hope) and 
will not be considered further here. The fun begins when new particles are 
produced in pairs in non-resonant processes, in particular those in which some 
particles go unobserved. Examples within the context of R-parity conserving 
supersymmetry might include processes like gg — > gg, gq — > gq, qq —>■ X2X1 
or qq — > 11. 

In this paper we present a new event variable, "mTGen" , which is designed 
to measure the mass scale(s) associated with any new particle(s) which might 
be pair produced at future colliders. On an event-by-event basis, mrcen sup- 
plies a lower bound for the mass of either of the two particles which were 
pair produced and whose decay products were observed, under the assump- 
tion that the event was indeed of that kind. The intention is that over many 
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events a histogram of m^ceii would reveal one or more edge structures whose 
upper endpoints would correspond to the masses of the particles which were 
being produced in large numbers. The rriTGen variable takes as input (1) 
the reconstructed momenta of the observed particles in each event, and (2) 
the masses of any unobserved particles which are taken to have been pro- 
duced in the decays of the primary particles. Though rriTGen can benefit from 
information regarding particle identification, this is not a requirement. 

It is not the intention of this paper to discuss how rriTGen performs when 
particle momenta are poorly measured, nor to discuss issues of finite precision 
or acceptance. Demonstration of the performance of m-TGen in environments 
representative of future colliders is the subject of current ongoing work. This 
paper seeks primarily (1) to provide a source of documentation for the defini- 
tion of the rriTGen variable, and (2) to document in the Appendix an analytic 
closed form approximation for which is valid for events with very little 
initial state radiation and which is needed to demonstrate that rriTGen can 
be calculated efficiently 

2 m TGe n 

As far as rriTGen is concerned, an event at a future collider is a set O contain- 
ing no observed Lorentz four-momenta: O = {of : i — 1, . . . , n }- Although 
quantum interference means that terms like "initial state" and "final state" 
cannot really be applied to the momenta in O in a well defined way, it is nev- 
ertheless common and expedient when analysing real events to treat some 
momenta as if they were the result of the decays from the pair of partons 
produced in the primary 2^2 process used in the matrix element, and to 
treat other momenta as if they were the result of initial state radiation (ISR). 
We will therefore divide the observed momenta into two non-overlapping sets 
F = {/f : i — 1, . . . ,Uf} (for momenta supposedly from the central 2^2 
pair production process) and G = {gf : i — 1, . . . , no} (for momenta suppos- 
edly from initial state radiation). In practice this might be done by assigning 
all four-momenta whose transverse momenta are greater than some threshold 
and whose rapidity is sufficiently central to F, while placing anything else in 
G. The value of rriTGen for each event will inevitably depend to some extent 
the exact nature of the cuts used, and potential biases would have to be in- 
vestigated in specific cases. Most event variables (e.g. thrust and sphericity) 
have similar second order dependence on cuts. The endpoint structures which 
are the eventual target of any rriTGen investigation, however, are expected to 
be particularly insensitive to these cuts. Individual events where misassign- 
ments are made will either be swept below the endpoint by the minimisation 
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procedure (when momenta are omitted from F), or will be smeared above 
the endpoint (when ISR with unusually large transverse momentum is added 
to F in error). 

In a real event, we do not know from which "side"0 of the event any 
particular observed particle has come. If we did know from which side each 
particle had come, we could use m T2 [U E] to place a lower bound on the 
mass of the two hypothesised outgoing primary particles. For a given event, 
even though we do not know the "correct" side assignments, there is nothing 
to prevent us trying all possible side assignments, evaluating for each 
of them, and then reporting the lowest value of itlt2 so obtained. This is in 
fact how rriTGen is defined. 



2.1 Definition of TUrcen 

m TGen is defined to be the smallest value of m T2 [U E] obtained over all 
possible partitions of momenta in F into two subsets a and (3 - each subset 
representing the decay products of a particular "side" of the event. Recall 
that my 2 is itself defined in terms of p^. and m a (respectively the transverse 
momentum and invariant mass of one side of the event), pj, and mp (re- 
spectively the transverse momentum and invariant mass of the other side of 
the event), and x (the mass of each of the unobserved particles which are 
supposed to have been produced on each side of the event) as follows: 
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where 
in which 



m 2 T (PT, Pr^ma, x)=m 2 a + x 2 + 2(E%E* 1 - p£ • p* 1 ) (2) 
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£t = \/(Pt) 2 + < and £^ = V(P?) 2 + X 2 (3) 



and likewise for a < — ► (3. With the above definition (in the case x = m x?)' 
m<T2 generates an event-by-event lower bound on the mass of the particle 



1 We will use the term "side" to refer to the division of the particles in F into two 
groups, depending on which of the two outgoing primary particles they descend from. An 
event, then, is an object with two sides, and possibly also some initial state radiation. The 
term "side" is not meant to suggest that the momenta of a particular side are in some 
way spatially correlated (e.g. in one hemisphere). Indeed, if the two primary partons were 
scalars produced at threshold, then the decay products of each "side" would be completely 
intermixed. 



3 



whose decay products made up either of the two sides of the event, under the 
assumption that the event represents pair production followed by decay to the 
visible particles and an unseen massive particle on each side. When evaluated 
at values of x m x° the above properties are retained approximately (see 
[HE]). There exist events which allow this lower bound to saturate, and so 
(in the absence of background) the upper endpoint of the m T2 distribution 
may be used to determine the mass of the particle being pair produced. 

We caution the reader to avoid the trap of mistakenly concluding that 
rriTGen, as defined above, is a function of purely transverse (and not also 
longitudinal) momenta of the visible particles in F. On the contrary, the 
definition above makes use of the z-momentum of every visible particle. Al- 
though this z-momentum plays no part in forming the transverse momentum, 
p T or pj,, of either side in any partition, it nonetheless can play a significant 
role in forming each side's invariant mass: m a or rap. (This caution is not 
unique to mTGen, but is equally relevant for any situation in which vtiti is 
used where a side consists of an "effective" particle composed of two or more 
real particles.) We note that it is possible to define a "Truly Transverse" 
form of m TGen , which we shall denote u m TTGen " by requiring, before evalu- 
ation begins, that each input four-momentum be individually longitudinally 
boosted to a frame in which its ^-component of momentum is zero. In effect 
this throws away each particle's z-momentum, does not change its transverse 
momentum, and reduces its energy so as to keep the particle's mass invari- 
ant. It is equivalent to evaluating mTGen in a "transverse" Minkowski space 
with (1+2) dimensions rather than the usual (1+3) dimensions. 

3 Discussion 

m<T2 has been used in the definition of rriTGem rather than a function of the 
invariant masses of the two sides of the event, as it is vital to account for the 
energy-momentum of unobserved particles. If an event has a non-zero total 
transverse momentum, then either there were unobserved particles in the 
event, or some momentum was mis- measured (or lost or created), or indeed 
a combination of both, ttiti takes these unobserved particles into account 
in hadron colliders. However, see the comment in the next section regarding 
use in lepton colliders. 

3.1 Regarding specialisations 

There is room for some slight specialisation of mTGen to the particular prob- 
lem in hand. mTGen need not always be calculated in exactly the same way. 
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For example, if it were desirable to suppose that neither side of an event 
could decay entirely into invisible particles, then one might wish to impose 
the additional constraint that F contain at least two momenta, and that one 
should only consider partitions of F into non-empty subsets. Events with 
fewer than two momenta in F could either be ignored or given rriTGen values 

of oE 

In an alternative specialisation, one could suggest constraining the parti- 
tions of F in such a way as to allow only those which meet certain require- 
ments in terms of conserved quantum numbers. For example, one might 
choose to reject assignments which make the absolute value of the charges 
on either side of the event greater than 1 or veto events where the total 
charge in F is more than some fixed valueJl One might try to restrict on the 
basis of lepton number, arguing from the standpoint of lepton universality 
(though this might be hard given the possibility of unobserved neutrinos). 

Finally, if mycen were to be used at a lepton collider where the momentum 
of the centre of mass was known to a reasonable precision, it would be sensible 
to replace m T2 (which is a variable designed for hadron colliders and so uses 
only transverse quantities in order to be insensitive to longitudinal boosts) 
with a variable analogous to itlt2 but designed to make use of z-momenta. 

3.2 Avoiding over-specialisation 

Though one can tailor rriTGen to the requirements or assumptions of a partic- 
ular investigation, it should be pointed out that the "philosophy" of rriTGen 
is to avoid such specialisation where unnecessary. 

Most events that are actually of the pair-produced type we are concen- 
trating on will have a large number of partitions. One and only one of those 
partitions is correct, and it will produce an m,T2 value that is appropriate 
(i.e. bounded above by the mass of the initially produced particle). All other 
(i.e. all wrong) partitions will result in two or more particles from different 
"sides" of the event being assigned to the same side. These wrong partitions 

2 When investigating particular classes of i?-parity violating supersymmetric models, an 
alternative specialisation might be appropriate. Here it might be desirable to suppose that 
none of the decay products on cither side of the event could go undetected (be invisible) . 
In this case one might consider evaluating not m,Ti in the definition of rriTGen but instead 
the larger of the two invariant masses of each side of the event. Note, however, that 
if this were done, one would have to satisfactorily address the question of where any 
observed missing transverse momentum had come from. Under the supposition at work 
here, large missing transverse momenta would indicate either large measurement error, or 
an event incompatible with the assumptions, and might suggest that the event should be 
disregarded or recalibrated. 

3 For example, at the LHC the charge of primary interaction should be —1, or +1. 
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thus tend to have have large ttit2 values, if only for the reason the union of 
particles from opposite sides of the event tends to yield four-momenta with 
very large invariant masses. Most wrong partitions therefore lead to val- 
ues larger than that generated by the "correct" partition, and so excluding 
large numbers of these "bad" partitions often has no effect on the eventual 
value of m TGen . 

By avoiding complex specialisations, one can make the variable to a large 
extent insensitive to a detector's ability (or inability) to distinguish particle 
types. This can make the variable useful for example during early running 
when the detectors' abilities to determine particle identification may not be 
well understood. 

3.3 Use as a cut variable 

Although m<T2 was originally proposed as a variable for measuring particle 
masses, it can also be used as a "cut variable" intended to separate certain 
new-physics signals from Standard Model backgrounds]^ This success is in 
part due an accidental conspiracy of three effects. (1) tends to small 
values for back-to-back QCD-like events. (2) tends to small values when 
the missing transverse momentum is small (which it is in much of QCD). 
(3) mrp2 tends to small values when the missing momentum is parallel to 
one of the visible particles fed to my 2 . This can easily happen in QCD 
when there are neutrinos in a jet, or a single jet is mismeasured through 
inadequate containment in a detector or passage through a crack region. By 
accident rather than by design, therefore, rriT2 tends to shift badly measured 
and Standard Model events away from large values, which are where the 
endpoints of the distributions containing new physics are expected to be 
found. 

T^TGen shares those features of m^, as it is always bounded above by m^2 
for the correct partition. m^Gen may therefore also be expected to find a role 
as a "cut variable", distinguishing events by their inherent mass scale, and 
able to focus on better-measured events. 

3.4 Comparison with other mass-scale variables: Mea? 

In the past it has been suggested that a good starting point for the de- 
termination of the mass scale of these new particles is the "effective mass" 
distribution j6j[7]. There are a number of slightly different definitions of Mes 

4 When used in this way the parameter \ which represents the mass of the stable invisi- 
ble particle is typically set to zero. This is of course correct to a very good approximation 
for the neutrinos which produce the missing momentum in Standard Model events. 
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and the phrase "mass scale" (a comprehensive list and comparisons between 
them may be found in [8]) but a typical definition of M Eff would be 



in which $t is the magnitude of the event's missing transverse momentum 
and where pru) is the magnitude of the transverse momentum of the i-th 
hardest jet or lepton in the event. 

All definitions of Mes are motivated by the fact that new TeV-scale mas- 
sive particles are likely to be produced near threshold, and so by attempting 
to sum up the visible energy in each event, one can hope to obtain an estimate 
of the energy required to form the two such particles. Broadly speaking, the 
peak in the Mes distribution is regarded as the mass-scale estimator. 

3.4.1 Problems with Mes 

Although the effective mass is a useful variable, and simple to compute, it 
has a few undesirable properties: 

Mes can be sensitive to the beam energy and the proton's parton dis- 
tribution functions (pdfs). This is primarily because the desired correlation 
between Mes and the mass scale relies on the assumption that the particles 
are produced near threshold. While it is true that the cross sections will usu- 
ally peak at threshold, they can have significant tails extending to \/I values 
considerably beyond the threshold value. As a result, the Mes distribution 
is broad, having a width similar to its mean. This smearing means that 
it is very hard to make precise statements about the mass scale from Mes 
alone. In contrast, the endpoints of m TGen are "precision measurements" 
in the sense that, up to the statistical error, the experimental resolution, 
and the systematics from mis- assignments between F and G, the endpoint 
of an mTGen edge has a direct interpretation as the mass of the underlying 
generated particles. 

The smearing due to pdfs described above means further that in the 
upper tails of the M E s distribution one has little chance of distinguishing the 
different contributions arising from distinct pair-production processes. There, 
the components of the Mes distribution are expected to blur together. With 
mTGen it is possible, though not guaranteed, that one will see multiple edges 
or changes in gradient in the upper parts of the spectrum allowing the masses 
of more than one high mass particle to be determined with precision. The 
extent to which this will work in practice is likely to depend on the extent to 
which the mj-Gen distribution is itself smeared by mis-assignments between 
F and G. Secondary edges at lower values of mrcen may, if not obscured by 




(4) 
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Standard Model backgrounds, also carry precision information on the masses 
of lighter particles. 

In fairness to M E s, one should point that mrcen has its own weaknesses. 
Most obviously, it is a number of orders of magnitude more costly^l to com- 
pute than Mes, for what could turn out to be little gain. Secondly, due to 
its m T2 dependence, m TGen has an explicit dependence on the hypothesised 
mass(es) x °f any invisible decay products. This makes it more complicated 
to handle, as it is should strictly be termed an event function (of x) rather 
than an event variable. It should be said in defence of mrcen, however, that 
it is good that this dependence on x is explicit and plays a precise role in 
offsetting the m^Gen endpoints to locations that have a clear physical inter- 
pretation. Though the M ES variable itself is simpler to compute, lacking x 
dependence, the need to consider x cannot be escaped when using M ES to 
draw conclusions about mass scales [8]. Thirdly and finally we note that 
mTGen, as canonically defined, may be more sensitive to high rapidity ISR 
than M ES (see section 13. 5[) . Having said this, the truly transverse form, 
m TTGen, should be much more tolerant of ISR while retaining the theoretical 
properties of mj- Gen J^l 

3.5 Example m^en distributions 

Simulations have been performed for several different supersymmetric parti- 
cle spectra, including the Snowmass points[H], for proton-proton collisions at 
LHC centre-of-mass energy of y/s =14 TeV. The HERWIGpH QU EE] Monte 
Carlo generator was used to produce inclusive unweighted supersymmetric 
particle pair production events. Final state particles (other than the invisible 
neutrinos and neutralinos) were then clustered into jets by the longitudinally 
invariant clustering algorithm for hadron-hadron collisions[13j used in the 
inclusive mode with R = 1.0 [13]. Those resultant jets which had both 
pseudo-rapidity (rj = — lntan6 l /2) satisfying \r]\ < 2 and transverse momen- 
tum greater than 10 GeV/c were used to calculate moreen and Mes- 

In figures [H [21 and [3] we show the distributions which would be obtained 
for several different spectra if it were possible to accurately assign all visible 
momenta to the correct category F or G (i.e. "interesting final state mo- 
menta" versus "initial state radiation"). The HERWIG initial state radiation 

5 Though see more positive view of this in section [3771 

6 Despite having ibttgen < titgen on an event-by-event basis, the position of the 
upper kinematic endpoint of the mxTGon distribution should be the same as the endpoint 
for mTGen (although with reduced statistics at the endpoint) provided that it is kinemat- 
ically possible for the decay products of each side to be produced with vanishing relative 
rapidity. 
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Figure 1: On the left hand side is a graphical representation of the susy 
mass spectrum of Snowmass point 4. The vertical positions of the particles 
indicate their masses. The horizontal positions of the centres of the bars in- 
dicate the relative LHC production cross-section (arbitrary units). The lines 
joining particles indicate decays with branching fractions in the following 
ranges: greater than 10 _1 solid; 10~ 2 — > 10" 1 dashed; 10~ 3 — > 10~ 2 dotted. 
The middle plot shows the distribution of our variable, rriTGen, with mrcen 
increasing vertically to ease comparison with the spectrum. The right hand 
plot shows the distribution of another variable, Meh/2, where Meat is defined 
in (TjJ. In both the m-TGen and the M Eff plots, the lighter shading shows the 
histograms with the number of events multiplied by a factor of twenty, so 
that the detail in the upper tail may be seen. 
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Figure 2: As for figure [TJ but for Snowmass point 5. 
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and underlying event have been switched off, and the parameter x which is 
required to calculate m TGen has been set to the mass of the lightest super- 
symmetric particle. The missing transverse momentum has been calculated 
from the negative vector sum of the momenta of the fiducial jets for reasons 
of computational efficiency as described in section 13.71 

It can be seen that the upper edge of the distributions gives a very good 
indication of the mass of the heaviest pair-produced sparticle. Distributions 
from a variety of different supersymmetric points show similar behaviour. 
This means that the position of the upper edge of the m-TCen distribution 
can be used to find out about the mass scale of any semi-invisibly decaying, 
heavy, pair-produced particles. 

Note that we have deliberately not used any information about the iden- 
tity of the observed particles, so we do not know from this plot alone whether 
the particles produced were squarks, gluinos, or indeed something completely 
different. But we do have a very good indication that there is a particle being 
pair-produced, and subsequently decaying to a mixture of visible and invisi- 
ble particles, and we have a good information about the mass scale at which 
this particle (or these particles) may be found. 

Furthermore, in all three change in slope can be observed at lower 

masses due to significant pair production of lower-mass particles (chargino 
and/or neutralino pairs for figure [Hand [3] or stop pairs for figure [2]). Therefore 
it is possible in principle to extract from this distribution information at 
several different mass scales. 

The plots also demonstrate some of the the undesirable properties of the 
variable Mes- There is, as has already been shown in [8], some correlation 
between M-^s and the mass scale of particles being produced. However Meb 
has a considerable tail at higher values caused by production of sparticles 
above threshold. 

The effects of Standard Model processes (and selection techniques re- 
quired to reduce them) are beyond the scope of this paper. More detailed 
studies using a complete set of Standard Model backgrounds and detailed 
detector simulation will be a important component of future work. How- 
ever we note that most SM events will have a small value of mTGen for the 
reasons discussed in section 13.31 Therefore, while these backgrounds might 
be expected to make it more difficult to extract information about lighter 
particles, we do not expect them to significantly affect the upper edge of 
the m-rGen spectrum which contains the information about the heavier (here 
squark and/or gluino) particle masses. 

The extent to which the approximations and assumption used in the 
calculation of mTGen can be justified is explored in figure HI If particles from 
initial state radiation and the underlying event are allowed to "pollute" the 
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Figure 3: As for figure [T] but for a (non-Snowmass) point with a heavier 
sparticle spectrum, defined by the mSUGRA parameters: {m = 1200 GeV, 
mi = 420 GeV, tan/3 = 10, m t = 174 GeV, \i < 0} and with a spectrum 
generated using Isajet[15] version 7.58. 
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events events events events 

(a) (b) (c) (d) 

Figure 4: mrcen distributions for the same spectrum as for figure El but 
with different assumptions, (a) The idealised described in text, (b) 

As for (a) but now with initial state radiation and the underlying event, 
and including particles from both categories F and G ("interesting" and 
"ISR/underlying event") to form jets, (c) As for (a) but with the parameter, 
X, (corresponding to the mass of the invisible particle) set to zero, (d) As 
for (a) but using the sum of the transverse momenta of the invisible particles 
for the missing transverse momentum. 
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final state (figure Hb) there is some smearing of the edgesQ The effect of 
adding the ISR and underlying event is to shift the apparent position of the 
m-TCen edge up by approximately 15-25% for our point (figure 0b). This value 
represents an upper limit to the uncertainty originating from these effects, 
which would be obtained in the unlikely case where the corrections from ISR 
and the underlying event were completely unknown. The extent to which 
this effect would be seen in experimental distributions will depend on the 
details of the event selection (for example on the rapidity cut on the jets). 
We note also that, as anticipated, if the same plots are generated for the 
"Truly Transverse" variant of the variable, rriTTGen, then the sensitivity of the 
endpoint to ISR is reduced (to around 10-15%) and there are proportionally 
fewer events at the endpoint. It remains to be seen whether the better 
strategy for the future will be to invest time in improving ISR rejection 
while focusing on mTGen > or to use variants like rriTTGen with less sensitivity 
to ISR but fewer statistics near the endpoint. 

If the invisible particle mass is unknown, or if m TGen is being used as a 
selection variable, then the distribution with the mass parameter x = is 
most appropriate (figure 0b). In that case the lower limit of the distribution 
(which cannot drop below x) is pulled down toward the origin. The change 
in the position of the upper edge is much smaller, as one would expect in 
the case where the total energy in the final state is dominated by visible 
particles. 

An efficient method of calculating mTGen which approximates the missing 
transverse momentum by the negative vector sum of the jet transverse mo- 
menta (see section \3.7h has been used for all these previous plots. This can 
be justified by its good agreement with the corresponding distribution ob- 
tained with the full numerical calculation of m TG en using the "true" missing 
momentum (figure [Hi). 

3.6 Cross section constraints 

mTGen relies purely on the kinematics of four-momentum conservation in 
each event. It makes no use of cross section information, which will therefore 
always remain a vital tool, orthogonal to mTGen; with which to constrain the 
overall mass scale. 

7 In the previous "unpolluted" plots, the assignment of final state particles to class F 
or G was determined from the kinematic history in the Monte Carlo event record. 
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3.7 Evaluating rriTGen 

Historically the main hurdle to the adoption of rriTGen has been the cost of 
evaluating it@ For each evaluation, m,T2 typically requires a numerical min- 
imisation to be performed which can take from a few thousandths to a few 
tenths of a second on most computers. The definition of mrcen requires this 
minimisation to be repeated up-to times: once for each partition of F. 
Since typical jet definitions can lead to high jet multiplicities in supersym- 
metric events, up to about 20 jets in some parts of parameter space, a naive 
implementation can take many seconds, or even hours to calculate rriTGen for 
a single event. 

For m TGcn to become usable, it is important to find less time-consuming 
ways in which the internal mj> 2 values can be calculated. Ideally, an ana- 
lytic or closed form expression for rriT2 is required to avoid time-consuming 
numerical minimisations. 

The authors were fortunate to be contacted in November 2006 by Ky- 
oungchul Kong and Konstantin Matchev (KKKM) [16]. KKKM informed 
the authors that they had derived an (undisclosed) analytic expression for 
rriT2 valid for the special case in which the missing transverse momentum 
of the event was entirely balanced by the transverse momentum of the two 
"key visible particles" input to rriT2- In other words, this special case cor- 
responds to having no net initial state radiation, or in the language of this 



paper 



= 0. It was not until June 2007 that the authors re- 

' T 

alised that this happens to be an interesting limit from the point of view of 
rriTGen- If all reconstructed momenta are thrown into F, then G is empty 
by construction, and the special case is satisfied. The more complicated the 
event is, the more grounds one has for doing this, as the less sure one can be 
as to the provenance of any individual particle. 

It is important to confirm the existence of the analytic expression claimed 
by KKKM and to perform timing tests using it in order to support the claim 
that rriTGen is now calculable in a reasonable time. The authors were not 
able to obtain the full expression for from KKKIvlJ^l It has therefore been 



8 The authors first proposed m^Gon for the analysis of the ATLAS Blind Data Challenge 
which concluded with the prize-giving at the ATLAS overview week in Prague, 2003, 
reported in New Scientist. See "Observations concerning the first ATLAS Blind Data 
Challenge" Barr, A J; Brochu, F M; Lester, C G; Palmer, M; Sabetfakhri, A; Aug 2003 
3 Though not releasing the full expression for m^, KKKM did release the answer for 

= one also has the masses 



g) If. 
= 1 Hi 



the special-special case where, in addition to 

of the two visible particles and the masses of the two hypothesised invisible particles all 
equal to zero. Though we did not make use of this information when deriving our own 
expression (it cannot be used for wtg™ as most partition create event "sides" which have 
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necessary to re-derive (what is hopefully) the same result independently in 
the Appendix. The authors understand that KKKM will release their own 
result in the near future [T7j. 

Using the analytic form of rriT2 derived in the Appendix, we find that even 
with a naive "try every partition of F" algorithm we can calculate mrcen for 
a 20-particle event in order one-second on a typical personal computer^ 
The computation time scales as 2 N where N is the number of particles, so 
a 10-particle event can be processed in one thousandth of that time. The 
authors find this is more than fast enough to make mrcen as usable as other 
standard event variables!"] 



4 Conclusion 

In conclusion, we believe that m^Gen could be an invaluable variable for 
physicists working at the LHC, and other future colliders. We hope that if 
its usefulness is validated in subsequent dedicated detector studies, it will 
become a standard tool in the Swiss army knife for new physics searches. 
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very large non-zero masses) we are grateful to have been able to use it as a check. 
10 This implementation is available from the authors on request. 

11 One would expect that, with more thought, it would be possible to create faster 
algorithms. Instead of trying all possible partitions when minimising over ttiti values, it 
might be possible to design an algorithm which flips momenta one-at-at-time from one 
side of the event to the other - hunting for the minimum. More work would be necessary 
to determine whether such algorithms could become stuck in local minima. 
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A Appendix 



Here we derive an expression for rriT2 for the special case in which the missing 
transverse momentum is entirely balanced by the two visible particles' trans- 
verse momentum - i. e. there must be no ISR. The particles themselves (and 
the hypothesised missing particles) are allowed to have arbitrary masses. 

Within this appendix we adopt the same conventions and definitions of 
[5] . In the language of that paper, the assumptions of our "special case" can 
be phrased as "S = a". We will begin, however, in the general case (E ^ a) 
and only introduce the simplification of the special case when we can no 
longer make progress without it. Beginning with the general case, then, we 
have the following notation: 



a" 


Lorentz 1+2 momentum of key visible particle 1 


(5) 




Lorentz 1+2 momentum of key visible particle 2 


(6) 




Lorentz 1+2 momentum of junk or ISR 


(7) 


p n 


Lorentz 1+2 momentum of invisible particle (mass x) produced with particle 1 


(8) 


<f 


Lorentz 1+2 momentum of invisible particle (mass x) produced with particle 2 


(9) 


A m 


Lorentz 1+2 momentum of unit a mass particle which is stationary in the lab fraliiS) 




real parameter (the reduced centre-of-mass energy from eq. 18 of [5]) 


(11) 



which are related by 

a" + ^ + gf + jf + (f = V^A». (12) 



Note that there is a potential ambiguity here between real momenta, mea- 
sured momenta, and hypothesised momenta. In this document, the quantities 
which are directly visible (a M , (3^ and g M ) are taken to be be real momenta, or 
equivalently to be measured quantities with zero measurement error. Con- 
versely g M and y/s are quantities which cannot be measured. In this 
case these symbols refer to the hypothesised neutralino momenta and/or hy- 
pothesised centre of mass energies that are used throughout the process of 
describing the event while attempting to calculate mx2- 
For simplicity, some derived quantities are also defined: 

Lorentz 1+2 momentum sum of the two key visible particles (13) 
Lorentz 1+2 momentum difference of the two key visible parti^M) 
Lorentz 1+2 momentum sum of everything seen in the detecto(15) 
Lorentz 1+2 momentum sum of the two invisible particles. (16) 

We already know that the particular momenta of p and q which need to 
be hypothesised to generate the value of fall into one of two categories. 



A*- 



IX -|- p' 



W = + gt* 
B» = if + 
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Either they are in a "balanced" configuration in which (a + p) 2 = ((3 + q) 2 or 
the value of is achieved for an "unbalanced" configuration in which this is 
not true. It is easy to determine whether a given set of momenta {a;^, g M } 
generate mx2 from a balanced or an unbalanced configuration, and also easy 
to determine what mx2 is for the unbalanced cases. We concentrate first, 
therefore, on the harder case of how to calculate the value of if it has 
already been determined that it occurs in a "balanced" configuration. 

A.l Balanced configurations 

In the "balanced configuration" , the value of m.T2 will be the minimum value 
of (a+p) 2 over all allowed values of provided that the following constraints 
are satisfied: 

X 2 (17) 
X 2 (18) 
W + q) 2 - (19) 

The first two constraints just put p^ and on mass shell. The final constraint 
is the one that makes the configuration mass balanced. When the constraint 
of equation (TIT?]) is satisfied we will refer to the both (a + p) 2 and (/3 + q) 2 
as M 2 . The approach we will take will be to assume a fixed value of yfs 
and then solve the above three equations for p^ and g M . We then explicitly 
minimise the resulting value of M by varying ^/s. The resulting minimum 
value of M is the value of mx2 we seek. 

From equation (1121) we can see that for fixed y/s, the value of B^ is fully 
determined: 

B» = v^A^ - W (20) 

and so the sum p^ + g M is fixed. We therefore choose to parametrise the 
three degrees of freedom which p^ and q^ have collectively by writing them 
in terms of an unknown Lorentz 1+2 vector 7^ as follows: 

p"=2 S " + 7" (21) 
<f = \b» - (22) 

Our stated intention of determining p M and g M for fixed \/s is therefore re- 
ally a requirement to determine the three components of 7 M , from the three 



P 



(a+pf 
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constraints in equations f|T7|) . (TTHT) and ffT9l . By substituting the two equa- 
tions above into equations f|T7|) . ffTSl) and ffT9l it is easy to show that the 
constraints on 7^ are equivalent to the following: 

7-5 = (23) 
7.a = ~A.(B + a) (24) 

7 2 = -i(5 2 -4 X 2 ). (25) 

The form of the above constraints motivates solving for 7 M as a linear combi- 
nation of the three linearly independent vectors er M and w^ 1 = t^ vi ' o v B T . 
Doing this, one finds two possible solutions: 



where 



fP' ± l r'\jH 2 + ^(B 2 -4x 2 ) (26) 



^ = (27) 



and 



ttr 

Each solution corresponds to a kinematic configuration which is a valid real- 
isation of original mass constraints, but the value of M will almost certainly 
be different in each case. Since our intention is to find m,T2, we will eventu- 
ally want to retain only the solution which gives the smaller value of M. As 
we do not yet know which solution that is, we retain both for the moment. 

We have now accomplished what we set out to achieve in step one. For 
fixed a/s we have defined B^ (with equation (1201) ). This value may be sub- 
stituted into equations (12T)]) . (|2T|) and (I2"H1) in order to find 7 M . In terms of 7 M 
we can then find the values of p M and (via equations (12T1) and ( |22l) ) which 
lead to the so called "balanced" kinematic structure in which both sides of 
the event have equal invariant mass M. All that now remains to do, is to 
minimise the value of M so-obtained over all allowed values of \fs. 

It is at this stage that we now move the "No ISR" special case that may 
be summarised as S M = cr 1 ' or equivalently as g M = 0. This change only 
affects terms with 5 M 's in them as these are the only quantities containing 
E M . Having made the substitution S' 1 — > cr M there is a substantial amount 
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of cancellation within the expressions in terms of which 7^ is defined (equa- 
tion f[2T)]) ). The net effect of this cancellation leaves M's dependence on y/s 
in the relatively simple form: 

M 2 = E + ± Ay / ( v / ^ - D) 2 - C 2 (29) 

for suitable values of the real quantities A, C, D, E and A which do not de- 
pend on ^Js. It is straightforward to show that the minimum of this function 
occurs when y/s takes the valued! 

= D + . ° (30) 

All that is needed to complete the evaluation of mx2 in this special case, 
then, is to determine the quantities A, C, D, and A. (E is not needed to 
calculate the value of ^fs which minimises M 2 Q) We can then evaluate ^fs 
in terms of these quantities, allowing in turn w^, 7^, p M and finally 
mT2 to be calculated. It may be shown that the values needed are as follows: 



A = 


i(A.(7) + i^[(a.A)-(A.a)(A.A)] ) 


(31) 


C = 




(32) 


D = 


A.a 


(33) 


and 




(34) 


A = 




(35) 


\g_\ 


where 




(36) 


J = 


Id 2 " (A.A) 2 
4|a| 2 


(37) 



and the quantity \a\ 2 is always evaluated in the lab frame 1^1 Note that the 
condition expressed in equation (I19p means that we can use a number of 

12 There is also a stationary point at y/s = D — C/\J 1 — -j? but it can be shown that 
this is always unphysical. 

13 Added in 2009: For completeness we note that E = \ 2 + \{ m2 a + - h® 2 + 
I I~tt^ [°' 2 (A.A) — (A.cr)((T.A)] . We note further that once the value of y/s obtained in 
([5(Tf is substituted into the expression for M in (j2"§)) one hnds that M^ in = E + AD + 
CVA 2 - A 2 . 

14 In other words \a\ 2 = (cr.A) 2 — a 2 . 
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different expressions to finally evaluate mx2- The simplest would be m\ 2 — 
(a + p) 2 or m\ 2 = {ft + q) 2 - However it is arguably nicer to preserve the 
explicit symmetry between the two sides of the event by instead evaluating 
the average of these two identical quantities. If this is done, one ends 
up with 



Note that it was the last of these forms which was used to generate the 
statement of equation (l2"9~j) . 

A. 2 Unbalanced solutions 

As discussed in [5] , the value of mx2 does not always arise from a configura- 
tion of hypothesised momenta in which both sides of the event have the same 
invariant mass. These unbalanced solutions arise if the momentum splitting 
which places one of the hypothesised neutralinos at the same transverse ve- 
locity, v t = pt/Ej*, as its visible "partner" (thereby minimising the invariant 
mass of that side of the event) causes the invariant mass of the other side 
of the event (which is then fixed by momentum conservation) to be even 
lower. This statement is generally true, and does not require the move to 
the £ = a special case considered in the section dealing with "balanced" 
solutions. Nevertheless, in order to write the full expression for the £ = a 
case we need to take these possibilities into account. 

A. 3 Putting all cases together 

We can now combine the two previous results into the following complete 
expression for m,T2 valid for events in which the missing transverse momen- 
tum exactly balances the transverse momentum of the two important visible 
particles (i.e. valid for the case £ = a also known as "no ISR"). 




(38) 



X 2 + ^« + ml) + -(a.B) + (A. 1 ). 



(39) 



(m a + xf 



iff 





(a + p) 2 or equivalently 

(ft + q) 2 or equivalently 

jd + \ (m 2 a + Wg) + \(v-B) + (A. 7) otherwise 



iff 



(40) 
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where g M = (\/x 2 + Iq| 2 > q) with q = — S —a and p' 1 = (\/X 2 + |p| 2 ? P 



with p = — S —f3 and in which p M and g M (which must not be confused 

with the entirely different quantities and g M just mentioned!) are defined 
by 

jf = ±B» + -f, (41) 
(f = \B»-Y, (42) 



-" H" ± tc^j H 2 + ^(B 2 -4 X 2 ) (43) 
(choosing the sign which leads to the smaller value of m^) in which 

w» = -== (44) 
v — w 

with 

= e ^ T a v B T (45) 

and 

iP = ~2 A -(^ + (J ) r( 5 2y M _ m^l . ( 46 ) 

and where we have taken in the "no ISR no junk" case 

B» = v^A^ - a" (47) 



having set 



= £> + , C (48) 



1 - ^ 



where 



A = I(A.(7) + i^[((7.A)-(A.(7)(A.A)], (49) 



c = V^l 2 + T' (50) 

D = A.a and (51) 

A = f^^Z (52 ) 
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in which 



(7 


2 _ 


(A. A) 2 


4! 


(7 


2 



j = m :,: ; (53) 



and where the quantity |a| = ^/ (o\A) 2 — <x 2 is the magnitude of the visible 
transverse momentum evaluated in the lab frame. 

A.4 Addendum in light of arXiv:071 1.4526 |IS1 

(This section was added in 2009. The purpose of the addition is to draw 
readers' attention to a nice result of relevance to this paper, published in a 
later paper by different authors [IE] • It is hoped that by placing the forward 
reference to [18] here readers may be more likely to use of the simplification 
that the result affords.) 

It was observed in [TH] that, after some simplification, equation (|40[) may 
be written in the following equivalent but shorter form (still valid only for 
the case S = cr also known as "no ISR") : 

{m a + x? iff {m a + x? > (/? + <?) 2 , 

{mp + xY iff (mp + x) 2 > (a + p) 2 , 



III: 



T2 



X 2 + A T + J(l + 2AT \_ mi )(A 2 T - ml m 2 ) otherwise 



(54) 

where the new quantity At is defined by At = \/ m 2 + a 2 J + f3 2 + a.(3 

and where all other quantities are defined as previously for ( 1401) . It will be 
noted that At is very closely related to the contransverse mass of [19]. We 
note that it is proved in [TBJ that the whole of (154"1) is invariant under simul- 
taneous equal magnitude but anti-parallel boost of a and (3 in the transverse 
Lorentz 1+2 plane. 
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